Study of Resonant Structures in a Deformed Mean-Field by the Contour Deformation 
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Solution of the momentum space Schrodinger equation in the case of deformed fields is being 
addressed. In particular it is shown that a complete set of single particle states which includes bound, 
resonant and complex continuum states may be obtained by the Contour Deformation Method. This 
generalized basis in the complex energy plane is known as a Berggren basis. The momentum space 
Schrodinger equation is an integral equation which is easily solved by matrix diagonalization routines 
even for the case of deformed fields. The method is demonstrated for axial symmetry and a fictitious 
"deformed 5 He", but may be extended to more general deformation and applied to truly deformed 
halo nuclei. 

PACS numbers: 21.60.Cs, 21.10.-k, 24.10.Cn, 24.30. Gd 



INTRODUCTION 



In nuclear physics, like in atomic physics, expansion of many-body wavefunctions on single particle bases, generated 

l practice. The newl; 

starts with the Berggren completeness [yj 0, 0, flil HH . The Gamow Shell Model has proven to be a promising 



by a suitable potential has been common practice. The newly developed Gamow Shell Model 



i particle bases, generated 



tool in assessing the structure of weakly bound and unbound nuclei. The Berggren completeness is a generalized 
completeness which treats bound and resonant states on equal footing. The completeness is given by a discrete sum 
over bound and resonant states and accomponied by an integral over a non-resonant continuum of scattering states 
with complex energy. A complete many-body Berggren basis may then be constructed from discretized single-particle 
Berggren orbitals where many body Slater determinants are constructed from the discrete bound, resonant and non- 
i resonant continuum orbitals. This is the philosophy of the Gamow Shell Model in full analogy with the standard Shell 
| Model using a harmonic oscillator basis. 

In this paper we address the problem of how to accurately calculate for resonances in deformed fields. Further, it 
+f ' is discussed how to obtain a complete set of states suitable for use in scattering/reaction problems, where spectral 
representations of Green's functions are of great interest, and in other many-body applications. The study of reso- 
nances in deformed fields has so far only rarely been considered, and for special cases. In Ref. [16j the solution of the 
angular momentum coupled Schrodinger equation was considered for a deformed Woods-Saxon. They diagonalized 
r* ■ the deformed Hamiltonian using a Berggren basis generated from the spherical Woods-Saxon potential in position 
. ' space, and compared with other methods such as an expansion in oscillator functions and a direct solution of the 
coupled equations. In Ref. energy levels and conditions for bound states to become resonances and resonances 
to become bound states were investigated for an axially deformed Woods-Saxon potential, by solving the the radial 
Schrodinger equation for coupled channels with outgoing asymptotics. However, the coupled channels method used 
in Ref. |l7j does not easily generalize to the non-resonant continuum. This implies that a complete Berggren basis 
in a deformed field is difficult to obtain, and all evaluated observables will become complex quantities unless the 
non-resonant continuum is taken properly into account. In Ref. p"s| a different approach was considered. Their aim 
was to propose a method to obtain scattering wave functions in the vicinity of a multi-channel resonance on the real 
axis, then calculate the phase shifts, and investigate whether a resonance condition is met. Further this method allows 
for evaluation of observables where the continuum is properly taken into account, and the observables become real 
quantities. 

In this paper we propose an alternative method, starting with the momentum space Schrodinger equation given in 
Eq. 0. In Ref. @ it was shown how a complete set of Berggren states may be obtained by an analytically continuation 
of the momentum space Schrodinger equation in the complex fc-plane by utilizing the Contour-Deformation-Method 
(CDM). By a suitable choice of deformed integration contour L + we demonstrated that all physical resonances con- 
verges extremely fast with respect to number of integration points. Further it was shown that for a particular type 
of contour L + , stable solutions of all physical scattering amplitudes may be obtained by a spectral representation of 
the Green's function. The main difference between a momentum space approach and a position space approach (see 
e.g. Ref. |T^). lies in their different discretization schemes. In momentum space, it is the Bessel completeness which 
is discretized, while in position space it is the completeness of the one-body problem (for example a Woods-Saxon 
completeness) which is discretized. The obvious advantage of the momentum space approach lies in its immediate sim- 



2 



plicity. Firstly, the boundary conditions are automatically built into the integral equations. Secondly, the discretized 
Schrodinger equation is a complex symmetric matrix which is easily diagonalized, and last but not least convergence 
is obtained by increasing the number of integration points. 

In sec.[n]the 1-dimensional angular momentum coupled integral equations in momentum space are derived. Sec. lIIII 
briefly discusses how the integral equations may be analytically continued in the complex fc-plane by CDM, and the 
relevant equations for numerical implementations are given. In Sec. lIVI a deformed field of Gaussian type is introduced, 
and a brief study of different axially symmetric deformations and multipoles is given. Sec. derives the multipole 
components of the Gaussian potential in momentum space. Sec. I VII gives results for single-particle resonances in the 
deformed Gaussian potential, and finally conclusions are given in Sec. IVIll 



II. MOMENTUM SPACE REPRESENTATION OF THE SCHRODINGER EQUATION. 

The momentum space single-particle Schodinger equation is given by, 

|^„(k) + J dk' F(k,k')Vn(k') = E n ip n (k). (1) 

Here the notation ip n (k) = (k|^„) and (k|V|k') = V(k, k') has been introduced. The potential in momentum space 
is thus a double Fourier-transform of the potential in coordinate space, i.e. 



V(k,k') = f-^ J dr J dx' e- ih T V{r,r')e ik ' T ' 



(2) 



Here it is assumed that the interaction potential does not contain any spin dependence. Instead of an differential 
equation in coordinate space (integro-differential equation for non-local potentials), the Schrodinger equation has 
become an integral equation in momentum space. This has many tractable features. Firstly, most realistic nucleon- 
nucleon interactions derived from field-theory are given explicitly in momentum space. Secondly, the boundary 
conditions imposed on the differential equation in coordinate space are automatically built into the integral equation. 
And last, but not least, integral equations are easy to numerically implement, and convergence is obtained by just 
increasing the number of integration points. Instead of solving the three-dimensional integral equation given in Eq. 
an infinite set of 1-dimensional equations can be obtained by invoking a partial wave expansion. To this end the wave 
function ip n (k) is expanded in a complete set of spherical harmonics, i.e. 

ip n (k) =^2^nlm(k)Yi m (k), ip n i m (k) = J dkYi m {k)^ n {k). (3) 

lm 

By inserting Eq.[3]in Eq. (JTJ, and projecting Yi m (k) from the left, the three-dimensional Schrodinger Eq. (JTJ is reduced 
to an infinite set of 1-dimensional angular momentum coupled integral equations, 

— k 2 - E nlm I ipnlm(k) = - ^ / dk'k' Vi m ,l'm'{.k,k')ipnl'm'{k'), (4) 
^ ' I'm' J 

where the angular momentum projected potential takes the form, 

= [dk I dk'Y,* m (k)V{ky)Y Vmt (k'). (5) 



Here dk = dflsinf? dip. In many cases the potential is given in position space, so it is convienient to establish the 
connection between Vi m ,l'm'(k, k') and Vim t i' m >(r, r'). Inserting position space completeness in Eq. © gives 

Vi m , Vm '{k,k') = J dr J dr' JdkJ dfc'^; i (fc)(k|r)(r|I/|r')(r'|k')y ;m (fc') = 

JdvJ dr' ||dfeF4(fc)<k|r)}(r|F|r')|y dk' Y lm {k'){v'\k') } . (6) 

Since the plane waves depend only on the absolute values of position and momentum, |k|, |r|, and the angle between 
them, 8kr, they may be expanded in terms of bipolar harmonics of zero rank |Tflj . i.e. 



e ik-r 



= ^Y J ^Ji(kr)(Y l {k)-Y l {f)) =J2(2l + l)i l ji(kr)Pi(cos9 kr ). (7) 



(=0 1=0 
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The addition theorem for spherical harmonics has been used in order to write the expansion in terms of Legendre 
polynomials. The spherical Bessel functions, ji(z), are given in terms of Bessel functions of the first kind with half 
integer orders [2flLl2lj. 

Inserting the plane- wave expansion into the brackets of Eq. © yields, 

/ 1 \ 3/2 

dfcr^(fc)(k|r> = i^—j 47r»- , j,(Ar)^(0. 
/ 1 \ 3/2 

dk'Y lm (k')(r'\k'} = (— J 47ri I jV(fcV)Fi/ m /(»»). 
The connection between the momentum- and position space angular momentum projected potentials is then given by, 
V lmtl , m ,{k,k') = -i 1 '- 1 drr 2 dr'r' 2 ji{kr)V lm , Vm -{r,r')3 V {k'r'), (8) 

I" J0 Jo 

which is known as a double Fourier-Bessel transform. The position space angular momentum projected potential is 
given by, 

VWm'(r,r') = J dr J dr' ^(f)^(r,r')^ m '(^)- (9) 

No assumptions of locality/non- locality or deformation of the interaction has so far been made, and the result in 
Eq. (JSJ is general. In position space the Schrodinger equation takes form of an integro-differential equation in case of 
a non-local interaction. In momentum space the Schrodinger equation is an ordinary integral equation of the Fredholm 
type, see Eq. This is a further advantage of the momentum space approach as compared to the standard position 
space approach. 



III. ANALYTIC CONTINUATION OF THE MOMENTUM SPACE SCHRODINGER EQUATION BY 

CDM. 



In Ref. [9j we discussed and outlined a method which analytically continues the momentum space Schrodinger 
equation through the unitarity cut onto the second Riemann sheet of the complex energy plane. The method is based 
on deforming the integration contour and is known as the contour deformation (distortion) method (CDM). As shown 
in Refs CDM allows for accurate calculation of a complete set of single-particle states, involving all kinds of 

poles of the scattering matrix. However, in Refs. only spherically symmetric fields were considered. Here we 

wish to generalize the method to deformed fields, and therefore we write down the relevant equations for the most 
general case. The rules for analytic continuation of integral equation with general integral kernels are not outlined 
here, since they are the same as for spherically symmetric fields. We refer the reader to Ref. |9( for further details on 
analytically continuation of integral equations and CDM. 

The 1-dimensional coupled integral equations given in Eq. Q are analytically continued from the physical to the 
non-physical energy sheet by distorting the integration contour. Choosing a suitable inversion symmetric contour L + , 
as discussed in Ref. |£j, we end up with the analytically continued coupled integral equations, 

(^-k 2 -E n i m \i/j n i m (k) = -^2 J dk'k' 2 Vi m M m ,{k,k')tp n i, m ,(k'). (10) 

^ ' I'm' ^ 

Here both k and k! arc defined on an inversion symmetric contour L + in the lower half complex fc-plane, resulting 
in a closed integral equation. The index n represents a bound or resonant state. The eigenfunctions constitute a 
complete bi-orthogonal set, normalized according to the Berggren metric |3. llflllllll2Lll3lll4Lll5l |. In solving Eq. (|1()|> 
numerically, we choose a set of N grid points in k— space by some quadrature rule, for example Gauss-Legendre. The 
integral is then discretized by J dk — > Yli=i w i- On the chosen grid Eq. (|10fl takes a complex symmetric form for 
bound, resonant and non-resonant continuum states n 

^-kflpnlmii) y/ W i W ] k ikjVl m ^ ml {ki,kj)^ n Vm' (j) = E n itl) n i m {i). (11) 

j I'm' 
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Changing from a continuous to a discrete plane- wave basis, it becomes transparent that the coordinate wave function 
is an expansion in a basis of spherical-Bessel functions 

[2 N 

(f>nim(r) = y ~ Vwihji(kjr)ip n i m (i), (12) 

i—l 

where ip n i m (i) are the expansion coefficients. Defining the functions 

fi(hr) = ^l^y/mhji(hr), (13) 
and using the discrete representation of the Dirac-delta function 

^-^tUt (14) 

we get the expansion 

N 

4>nlm(r) =^2lpnlm(i)fl(kir), (15) 
t=l 

where it is easily seen that the functions fi(kir) are orthogonal for different fc, and normalized 

drr 2 f l {k i r)Mk j r) = 5 kiikj , (16) 



$ki kj being the Kronecker delta. The complete and discrete set of single-particle orbits defined by this contour will 
then include the pole states, i.e., anti-bound, bound and resonant states, and the discretized complex continuum 
states defined at each point on the contour. 



IV. DEFORMED FIELD OF GAUSSIAN TYPE. 



We consider an axially symmetric deformed Gaussian potential with no spin and tensor components. In polar 
coordinates it is given as 

V(r,6) = Vb exp (-r 2 (a cos 2 6> + /3 sin 2 (9)) , (17) 

or in Cartesian coordinates, 

V{x, y, z) = V exp (-f3{x 2 + y 2 )) exp(-az 2 ), (18) 

here Vb is the strength of the potential and a and (3 are shape parameters. Here z is the symmetry axis, and the 
potential is reflection symmetric in the x, y-plane. In the case a = (3 the potential is just a spherical Gaussian 
potential. In the case a > [3 the potential field is contracted along the z-axis, and defines an oblate shape. In the 
case a < [3 the potential field is stretched out along the z-axis, and defines a prolate shape. Defining a deformation 
parameter 8 by 

<5=l-~ (19) 

P 

Eq. 117|) can be written in the form, 

V(r, 9; (3, 6) = V exp(-/3r 2 ) exp(/3<5r 2 cos 2 9) = V(r; f3)D(r, 9; f3, S). (20) 

Here V(r; (3) is a spherically symmetric formfactor and D(r,8; (3,8) a deformation formfactor, D = 1 for 6 = i.e. 
a = (3. We require that the volume of the central potential, with the shape parameter «o = a = (3, is equal to 
the volume of the axially symmetric deformed ellipsoidal potential. This implies that the shape parameters of the 
non-central and central Gaussian potential satisfy the following relation, 



a/3 2 = al 



(21) 
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and the deformation parameter S may be expressed in terms of ao and by 

Fig- d shows plots of the isocurves V(r,8) = 0.5 in the x, z-plane for the deformation parameters 6 = ±0.5. With 
potential parameters ao = 1 and Vq = 1 for the spherically symmetric potential, S — 0.5 gives the shape parameters 
a = 2~ 2 / 3 and j3 = 2 1 / 3 for the deformed potential, and 5 = —0.5 gives the parameters a — (3/2) 2 / 3 and (3 = (2/3) 1//3 , 
respectively. It is seen that S — 0.5 corresponds to an prolate shape, for S = —0.5 the potential takes a oblate shape, 
the symmetry axis being the vertical z-axis. In order to assess the shape structure in more detail, it is instructive 
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FIG. 1: Plot of the ellipsoidal isocurves V(r, 9) = 0.5 of the deformed potential for deformation S = ±0.5 with potential strength 
Vo = 1 in the x, z-plane. For the spherically symmetric potential a shape parameter ao = 1 was chosen. 

to study the multipole components of the potential. An axially symmetric potential may be expanded in terms of 
Legendre polynomials, i.e. 

V(r,9) =J2 V x (r)P x (cos 6). (23) 

A 

Using the orthonormality properties of the Legendre polynomials, the multipole components are given by the integrals 

Vx(r) = (2A + 1 V q exp(-/3r 2 ) / drj exp(/3«5r 2 7 7 2 )P A (r;), (24) 



where 77 = cos((9). Here it is explicitly seen for our reflection symmetric potential, that only even multipoles give 
non-vanishing contributions, since the Legendre polynomials have the property 

Px(-v) = (-l) X Px(v), 

and the potential is an even function in 77. 

The monopole part of the deformed Gaussian potential may be calculated analytically, 

V x =o(r) = ^V exp(-(3r 2 ) J dr) exp(/MrV) = ^V exp(-/3r 2 ) A)(r), (25) 

where 

D (r) = ^=7(1/2, r), r = -pSr 2 . (26) 
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FIG. 2: Plot of A = 0, 2, 4 multipoles of the Gaussian potential with deformation parameter 8 = 0.5 and ao = 1. 
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FIG. 3: Plot of A = 0, 2, 4 multipoles of the Gaussian potential with deformation parameter 8 — —0.5 and ao = 1. 



Here 7(1/2, r) is the incomplete gamma function, see e.g. Ref. [23. Figs. |3 and |3| give plots of the A = 0,2,4 
multipoles of the Gaussian potential with deformation parameters 5 — ±0.5 and the potential parameters ao = 1 and 
Vo = 1. It is seen that the radial monopole distributions are more or less identical for S — 0.5 and S — —0.5. Further 
it is seen that the deformed Gaussian potential is nearly a pure quadrupole deformation, since the A = 4 multipole is 
almost vanishing in both cases. This may be understood from considering the exponent of the deformed formfactor 
in Eq. (J2DJ, which can be rewritten in terms of the Y 2 o(f) spherical harmonic. 
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V. MULTIPOLE COMPONENTS IN MOMENTUM SPACE. 



Having discussed the shape and multipoles of the deformed Gaussian potential, we now turn to the actual so- 
lution of the Schrodinger equation for this potential. We wish to solve the partial wave decomposed momentum 
space Schrodinger equation given in Eq. QJ, an d therefore need the Gaussian deformed potential in a partial wave 
decomposed, momentum representation. The Fourier transformation of the deformed Gaussian potential in Eq. 118(1 
is, 

V(q x , q y , q z ) = ° J dx dy dz exp (i(q x x + q y y + q z zj) exp (-(3(x 2 + y 2 ) - az 2 ) 

= 8,3/^1/2 6XP (-^(92 + «2>) exp ("4^) < W 
where qi = ki — k[i i = x, y, z. In terms of spherical momentum space coordinates k, 9, tp the potential takes the form, 

V(k,k') = exp | - ^(fc 2 sm 2 6> + fc /2 sm 2 6>') - — (fccosfl - fc'cosg') 1 

Stt 3 / 2 ^ 1 / 2 V 4 P 4a 

x exp ( -^-kk' sin 9 sin 9' cos (ip - tp') } . (28) 

Due to axial symmetry the dependence of the potential on the azimuthal angles ip, tp' is only on the difference 
lu = tp — tp' . The potential may therefore be expanded in a complete set of harmonics, i.e. 

oo 

V(k,k')= ]T ^(M')exp(i/iw), (29) 

here k = (k,9). The harmonics exp(i/xw) obey the orthogonality relation 

dus exp(— ifju) exp(i//a;) = (30) 
the /x'th harmonic of the potential is therefore given by the integral 

V^{k,k') = ^-j du exp(-^)^(k,k'). (31) 
From Eq. I|28[) it is seen, that for the integral over u>, we have to consider the following integral, 



i r i r 

Iu(y) — — / did exp(— ifju) exp(ycosuj) = — du cos(/Uw) exp (y cos w), 

271 " J -IT ' 1" JO 



(32) 



where we have introduced the variable 



y = — kk! sin 9 sin 9'. 

Ip 

The integral in Eq. (|32|l is just the definition of the modified Bessel function of the l'st kind (see e.g. HJ)- The /x'th 
harmonic of the potential is thus of analytic form, and given by, 

v fl (k,k')= hi 



87r 3 / 2 /3a 1 /2 

exp f — ^-(fcsinfl- fc'sinfl') 2 - — (kcos9- k' cos9'f ) exp(-y)IJy). (33) 
\ 4/J 4a / 

Inserting the expansion of the potential given in Eq. (|29|l into the momentum space Schrodinger Eq. and projecting 
the equation on the harmonics exp(ifjio), the three-dimensional integral equation has been reduced to an infinite set 
of two-dimensional integral equations. The /z'th integral equation is easily solved as a matrix diagonalization problem 
with dimension N r x No. Here N r is the number of integration points for the radial integral and Ng the number of 
integration points for the angle integral. However, the Schrodinger equation can be further reduced to a coupled set 
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of one-dimensional integral equations by projecting on spherical harmonics (see Eq. (@J). 
projected potential in Eq. 1(28(1 then takes the form, 

V lm>lw (k,k') = JdkJ dfc'l^(fc) j jr ^(MOexp(w)j^'m'(fc') 



The angular momentum 



2tt / d9 sin 9 I d6' sin6'Pi m {cos6)V m (k,k')P Vm {cas6')5 n 
Jo Jo 



where Pi m {x) are the normalized associated Legendre polynomial, 

- , , (2l + l(l~ m y.} 1/2 



(34) 



(35) 



In our calculations we start with the angular momentum projected potential given in Eq. 134(1 . In numerical 
calculations the multipole expansion in angular momentum I has to be truncated at some i max and the Hamiltonian 
matrix to be diagonalized is, 



H' 



H m *{hM) 



i/ m7r (7i, / max ) 



Timix (J 1 \ timTT// 7 

11 \ 'maxj «\) • ■ • 11 v max? l r\ 



(36) 



Here m and 7r are the angular momentum projection and parity, respectively, which are good quantum numbers in 
case of axial symmetry. The allowed values of / are even and odd for positive and negative parity states, respectively. 
Each submatrix H mn (l, I') in Eq. ((36(1 has matrix elements given by, 



Tjrn-K 
U i,j 



\JWiWjk{kjVi m ii m (ki^ kj). 



(37) 



The rank of each submatrix H m7T (l,V) are determined by the total number of integration points used in the dis- 
cretization of the integration contour L + in the coupled momentum space Schrodinger equation, given in Eq. 1(10(1 . 
The results reported in this work, used the same number of integration points for each coupled equation, so the total 
rank of the matrix to be diagonalized is N x AT; , where N is the total number of integration points and Ni the total 
number of angular momentum coupled integral equations given in Eq. ((10(1 . Diagonalizing the complex symmetric 
matrix ((36(1 . we obtain a complete set of states within the chosen discretization space. The basis may be utilized in 
different spectral represenations used in scattering theory or in Gamow-Shcll-Model calculations involving deformed 
fields. 

In all calculations reported below, we used a discretized contour L + defined by a rotation 9 and a translation C 
in the complex fc-plane, see Fig. 0] In Ref. it was shown that this type of contour allows for a stable numerical 
solutions of the scattering amplitude (or t— matrix), by using a spectral representation of the Green's function. In 
physical scattering, the energy is given along the real axis. By defining a basis with continuum energies given along 
the contour depicted in Fig. the problem with poles of the Green's function are eliminated, see Ref. |9j for more 
details. The total number of integration points is given by N — Nr + Nt, where Nr is the number of points along 
the rotated line defined by 9 and Nt is the number of points along the translated line defined by C. In numerical 
calculations fc max (see Fig. should be chosen large enough, so that the calculated wave functions and energies do 
not change with increasing fc max - In our calculations we used fc max = 6fm~ 1 . 



VI. FORMATION OF SINGLE-PARTICLE RESONANCES IN A DEFORMED GAUSSIAN 

POTENTIAL. 

As a model study we consider the Gaussian potential given in Eq. I(17|) . which in the spherically symmetric case 
reproduces the J x = 3/2^ resonance in 5 He. The J 71 ' = 3/2^ resonance, to be associated with the single- particle orbit 
P3/2> is experimentally known to have a width of T ss 0.60 MeV. 

In our calculations we used the following parameters for the spherically symmetric Gaussian given in Eq. 117(1 . 

V = -53.5MeV, a Q = 0.188fnT 2 . (38) 

As a check of our momentum space approach, we compared our results for the spherical limit with the results 
obtained with the computer program GAMOW |2^|. For the spherical Gaussian potential, GAMOW gives a bound 
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FIG. 4: Sketch of the contour L + in the complex fc-plane. The contour is specified by a rotation angle and a translation C 
in the fourth quadrant of the complex fc-plane. 



state for the V = 0+ channel with energy E = — 14.9044MeV, and a resonance for the V = 1~ channel with energy 
E = 0.7268 — 0.3096iMeV. Here the nucleon spin s = 1/2 is neglected since the energy levels are degenerate for 
j = I ± 1/2, which follows from the spin independence of the Gaussian potential. Our momentum space calculations 
were able to exactly reproduce these results with a total number of discretization points N = 25. This comparison of 
two different methods, provided us with a check of our codes and our derivation of the momentum space equations. 

In order to obtain converged results for the deformed case, we investigate the convergence with respect to total 
number angular momentum coupled equations, see Eq. (|36[) . and with respect to the total number of integration 
points used in discretization of each coupled integral equation given in Ea. (|10f> . First, convergence with respect to 
total number angular momentum coupled equations in Eq. (|36|1 is considered. We fixed the number of integration 
points to N — Nji + = 20 + 30 = 50, this is large enough to ensure convergence with respect to number of 
discretization points. For bound states we used a real contour L + , i.e. 8 — and the real k— axis was discretized 
with 50 points. In the case of resonant states we used a complex contour L + defined by a rotation 9 = ir/4 and a 
translation C = sin(7r/4) x 0.4fm _1 ~ — 0.29fm in the complex fc-plane (see Fig.QJ. 

Table 0] gives the convergence of the — + ground state energy, for deformation parameters 5 = ±0.9. For 
S = —0.9, convergence is quickly reached, with l max — 4. For 6 — 0.9, convergence is considerably slower. It is seen 
that the deformation S = 0.9 affects the bound state most, and the ground state becomes less bound E = — 12.1McV, 
for the prolate deformation. On the other hand, the oblate deformation S — —0.9 has little effect on the ground 
state energy E = — 14.7MeV. This may be understood by considering the monopole term of the potential, which is 





S = 0.9 


6= - 


0.9 


^max 


Re[E] 


Im[E] 


Re[E] 


Im[E] 





-10.5843 


0. 


-14.4816 


0. 


2 


-11.9041 


0. 


-14.6533 


0. 


4 


-12.0741 


0. 


-14.6551 


0. 


6 


-12.0953 


0. 


-14.6551 


0. 


8 


-12.0979 


0. 


-14.6551 


0. 


10 


-12.0983 


0. 


-14.6551 


0. 



TABLE I: Convergence of groundstate, m n = + , for deformation parameters 5 — ±0.9 as the number of partial waves increases. 
In the spherically symmetric case 8 — the l w = + Gaussian potential supports a bound state at energy E — — 14.9044MeV. 



the main component in the multipole expansion in Eq. H'-it . In Fig. [S] a plot of the monopole part of the Gaussian 
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r [fm] 

FIG. 5: Plot of the monopole part of the Gaussian potential with deformation parameters 5 = and S = ±0.9. 

potential with deformation parameters S — ±0.9 is given, together with a plot of the spherically symmetric potential. 
It is seen that the monopole term for the S — —0.9 potential is more or less identical to the spherically symmetric 
potential (slightly less attractive), on the other hand the monopole term for the S = 0.9 potential is less attractive for 
r < 4fm, but more attractive at large distances r > 4fm. From this one may conclude that the ground state of the 
5 = —0.9 potential will be more bound than for the S — 0.9 potential, since the ground state is deeply bound and the 
wave function will be mainly located in the interior part of the potential. 

The resonant orbit V = l~ in the spherically symmetric potential is split into two non-degenerate orbits ( m* = 0~ 
and rri K = 1~) in the case of an axially symmetric deformation. This is a characteristic of axially deformed potentials. 
Table [H] gives the convergence of the = 0~ and m" — 1" excited negative parity states in the Gaussian potential, 
for deformation parameters 5 = ±0.5. In all cases a satisfactory convergence is obtained with i max = 5. For the 
states with vanishing angular momentum projection along the z-axis ( m = ), it is seen that for 5 = 0.5 ( prolate 
deformation ) the V = 1~ state has become a bound state with energy E = — 0.680MeV. For zero angular momentum 
projection, the particle moves in an orbit along the z-axis. So in the case of a prolate deformation, where the field 
is stretched out along the z-axis, a particle moving in this orbit will "feel" the field more strongly than compared 
with the spherically symmetric field, and it will become more bound. This explains also why the particle with m = 
becomes more unbound in the case of the oblate deformation S = —0.5, see Columns 6 and 7 of Table ITU 
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m w 
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^max 


Re[E] 


Im[E] 
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Re[E] 


Im[E] 


Re[E] 


Im[E] 
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-0.5282 


0. 


1.4865 


-1.0177 


1.5402 


-1.0701 


0.3815 


-0.1139 


3 


-0.6772 


0. 


1.4419 


-0.9631 


1.5170 


-1.0404 


0.3602 


-0.1042 


5 


-0.6802 


0. 


1.4410 


-0.9621 


1.5168 


-1.0402 


0.3601 


-0.1041 


7 


-0.6803 


0. 


1.4410 


-0.9620 


1.5168 


-1.0402 


0.3601 


-0.1041 


9 


-0.6803 


0. 


1.4410 


-0.9620 


1.5168 


-1.0402 


0.3601 


-0.1041 



TABLE II: Convergence of the m n = and m w = 1 energies, for deformation parameters 8 = ±0.5 with increasing number 
of partial waves. In the spherically symmetric case (5 = 0) the J lr = 1~ Gaussian potential supports a resonance state at energy 
E = 0.7268 - 0.3096iMeV. 



For the m = 1 case the opposite takes place. In the case of S = 0.5 the particle becomes more unbound, while for 
S = —0.5 the particle becomes more bound. By considering the dipole (I = 1) term of the wave function, the particle 
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2.E-06 
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7 
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0. 


3.E-07 


-2.E-07 


2.E-08 


-l.E-08 
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-4.E-09 


9 


8.E-09 


0. 


4.E-09 
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9.E-11 


-7.E-11 


4.E-11 
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TABLE III: Convergence of the m 7T — and rri* = 1 squared amplitudes of the wave functions for each partial wave I, for 
deformation parameters 8 = ±0.5. 



moves in an orbit making 7r/4 degrees with the z-axis. From this it may be understood that the particle gains more 
binding in the case of an oblate deformation S = —0.5 and becomes more un-physical in the opposite case S = 0.5 ( 
see columns 4,5,8 and 9 of table IXTfl . 

In table ITTT1 the squared amplitudes of the wave functions are given for each partial wave I. It is seen that in all 
cases that the squared amplitudes for the I = 1 component of the total wave function, is nearly equal to the norm of 
the total wave function, while all other partial wave amplitudes are vanishing small. In this sense one may say that 
the orbital angular momentum is approximately a "good" quantum number. 

Having investigated the convergence with respect to number of angular momentum coupled equations, we now 
consider convergence with respect to number of discretization points along the contour L + , for the negative parity 
states. For the energies in Tables [H] and IIIII we reached satisfactory convergence with l max — 7. In considering 
convergence with respect to integration points, we then fix the maximum number of coupled equation to Z max = 7. 
Table iTVl reports the convergence of the odd parity energies in the deformed Gaussian potential given in Eq. 
with deformation parameters S = ±0.5. It is seen that one obtains convergence with a total number of integration 
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-0.6777 


0.0000 


1.4401 
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-1.0448 


0.3577 


-0.1091 


10 
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-0.6803 


0.0000 


1.4411 


-0.9621 


1.5170 


-1.0403 


0.3601 


-0.1043 


10 
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1.4410 


-0.9620 


1.5168 


-1.0402 


0.3601 


-0.1041 


10 


25 


-0.6803 


0.0000 


1.4410 


-0.9620 


1.5168 


-1.0402 


0.3601 


-0.1041 


15 


25 


-0.6803 


0.0000 


1.4410 


-0.9620 


1.5168 


-1.0402 


0.3601 


-0.1041 


15 
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-0.6803 


0.0000 


1.4410 


-0.9620 


1.5168 


-1.0402 


0.3601 


-0.1041 


20 


30 


-0.6803 


0.0000 


1.4410 


-0.9620 


1.5168 


-1.0402 


0.3601 


-0.1041 



TABLE IV: Convergence of the rri" = and m 7T = 1 energies with increasing number of discretization points along the 
contour L + . Here deformation parameters 8 = ±0.5 were used, and i max = 7. 

points given by N = Nr ± Nt = 10 ± 20 = 30. Note also that with N = 15 points, we have satisfactory results. For 
N = 30 and Ni = 4 (only four coupled equations for Z max = 7 due to conservation of parity) the total dimension of 
the matrix in Eq. (|36|l is dim = 120, which is diagonalized extremely fast with any diagonalization routine suitable 
for complex symmetric matrices. 

Having determined convergence properties of the odd and even parity states in the deformed Gaussian potential, 
we now study how the different states behave over a large range of deformations. In figure 10 a plot of the bound 
state energy of the m 11 = + state is given for the deformation parameter S taking values between — 0.9and 0.9. It is 
seen that the position of the bound state varies much more strongly for a prolate deformation (5 > 0), than for an 
oblate deformation. 

Fig. shows a plot of the real (solid lines) and imaginary part (dashed lines) of the rri" — 0~ and m* = 1~ states 
for the deformation parameter S taking values between — 0.9and 0.9. 

The value for 5 in which the m" = _ resonant state becomes a bound state is given when the real energy trajectory 
meets the imaginary energy trajectory for Re[i?] < 0, i.e. Im^] = Re[_E]. Here the splitting of the = 1~ resonant 
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FIG. 6: Bound state trajectory for the ml' = + state in the deformed Gaussian potential. Energy is plotted as function of 
deformation parameter S. 



level with respect to the angular momentum projection m is clearly seen. It is also seen that the energy of the 
m* = 0~ and the m 71 " = Instate behave in opposite manner for 5 > and for 6 < 0. The m r! = 0~ resonance for S = 
becomes a bound state for S > 0.3. For S < the to 71 " = 0~ resonance moves further down in the lower half complex 
/c-planc. On the other hand, it is seen that the rri K — resonance state does not become a bound state for the values 
of S considered, 5 e (—0.9, 0.9). For S e (—0.9, 0.2) the m 7 ' = 1~ resonance energy displays a weak variation from the 
6 = energy, and slowly moves towards the scattering threshold E = 0, for 5 — > —0.9. On the other hand, as 5 — > 0.9 
the imaginary part of the energy dives into the lower half complex energy plane, and the resonance state becomes 
strongly unphysical. 
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FIG. 7: Real (solid lines) and imaginary parts (dashed lines) of the m n = and the m Tr = 1 state energies in the deformed 
Gaussian potential as the deformation parameter S is varied between —0.9 and 0.9. 



13 



VII. CONCLUSION. 



In this work, we have demonstrated that the Contour Deformation Method may be generalized to the case of 
deformed potentials. CDM is applied to the momentum space Schrodingcr equation, allowing for stable and converged 
solutions of physical resonances. In addition a complete set of Berggren states is obtained, which may be used in 
the construction of a many-body Berggren basis. The most obvious advantage of this momentum space approach, as 
compared to its position space analog, is that the boundary conditions of all kinds of states, i.e. bound, resonant and 
continuum, are automatically taken care of since we are dealing with integral equations instead of integro-differcntial 
equations. The method is demonstrated for axial symmetry and a fictitious "deformed 5 He", but may be extended to 
more general deformation and applied to truly deformed halo nuclei. Such applications as to 8 Li, 9 Be and 11 Be will 
be discussed elsewhere. 
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